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ABSTRACT 


An expansion of the planetary disturbing force is given in a 
form convenient for a method of successive approximations in the 
application of the method of rectangular coordinates. The diffi- 
culties associated with slow convergence along powers of the plan- 
etary masses are absorbed in the iteration process. 
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COMPUTATION OF GENERAL 
PLANETARY PERTURBATIONS, PART III. 
AN EXPANSION OF THE DISTURBING FORCE* 

by 

Lloyd Carpenter 
Goddard Space Flight Center 


INTRODUCTION 

The present level of interest in general planetary perturbations emphasizes the importance 
of any developments which have strong computational advantages. One of the main difficulties in 
the computation of perturbation theories is the slow convergence along powers of the planetary 
masses. This slow convergence comes from small divisors associated with near resonance or 
possible close approaches of the planets or relatively large factors from the eccentricities. In 
such cases the convergence can become arbitrarily slow to the point where a perturbation theory 
is not appropriate. However, the convergence along powers of the masses is slow in many cases 
where the total perturbations are not too large. These ideas were borne out by the work of Hill 
(Reference la or lb) where the higher order terms were much larger than expected. In the 
method of rectangular coordinates (Reference 2) this difficulty can be reduced by modifying the 
usual concepts of order and replacing the expansion in powers of the disturbing masses by an 
expansion in powers of the perturbations themselves. 

In terms of the expansion this means that in the disturbing force the orders of terms are 
based on the powers of the perturbations and the mass factor, but after the integration all terms 
are treated alike or on the basis of the actual sizes of the coefficients. The equations are inte- 
grated by successive approximations since the unknown perturbations appear in the integrands. 

In the iteration process the start may be based on elliptic motion or any approximate or previous 
general theory. At any stage the coefficients may be modified to improve the convergence of suc- 
cessive approximations if necessary. 

The results may be compared in part with Musen (Reference 3). 


*In Part I (Reference 4) of this report, a description of a program for computing Hansen’s Planetary Perturbations was given. Part II 
(Reference 5) contains a comparison of Hansen’s components with the rectangular components in Musen s method. 
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THE BASIC EXPANSION 


Take the equation of motion in the form 


d 2 r 
dt 2 




i — 1 



r . 



The notations are 

r the position vector of the disturbed planet with respect to the sun 
r. the position vector of the i th disturbing planet with respect to the sun 

Jo. the position vector of the i th disturbing planet with respect to the disturbed planet (Jo. = r. - r ) 

M the mass of the sun 

m the mass of the planet 

m. the mass of the i th disturbing planet 

ij ? = k 2 (M + m) where k is the Gaussian constant 

t the time. 


It is understood that there are N disturbing planets. The equation of motion is decomposed into 
the equation of undisturbed motion 


plus the variational equation 



- M 


r 

r 


o 

3 

0 


with 


and 



F 



r = r n + s , 


where r 0 is the elliptic motion and s is the perturbation. Each of the r. and p i are separated in 
the same way. To obtain a useful form for computing general perturbations, the terms r/r 3 ? pVp? 
and r./r 3 are expanded in powers of the perturbations. Each of these terms has the same form. To 
simplify the development the standard binomial expansion is used in place of the elegant symbolic 
operators used by Musen (References 2 and 3). 
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Take the general vector q in the form 


q n + e 


Then 

q 3 =q 3 0 [l + S]3/2, 


where 


8 



e = |e | and q = |q| 


When S is small the binomial expansion 


1 

q 


3 






or 


1 = J. _3_j>_ 15 35 j >3 315 § 4 

q 3 q 3 2 q 3 8 q 3 16 q 3 128 q 3 


is useful. Substituting the expression for 8 


q q 0 ft q ° ‘ e 
q3 q o \ q o q o 


, q o ■ e 3 e 2 - 15 ( q o ' ?) 2 _ 

3 7- ? + v q °-7- — q ° 

L q o q o q o 


x j3 6 2 7 _ 15 ^o'^) 2 ? 15 q 0*^ 2 r 35 (^o*^) 3 


L 2 q s o 2 < 


2 q S 


q 0 +■ 


2 qj 


^0 


15 e 4 105 ( q o' e ) 2 fi2 315 ( q o ' e ) 4 

8 "a 2 + ^ 8 

q o q o q 0 J 


35 (%• O _ _ 15 

e ~ 2 


q ° + T 




The terms have been separated according to the various powers of t . 



It is convenient to separate the perturbation 7 into classes of terms. For example in the 
next section these classes will represent terms factored by various powers of the time, and in the 
following section they will represent contributions from various disturbing planets. Put 


and 




i , j . k 


Each component 7 k is treated as if it were of the first order (the same order as 7 ). Now 


7 - q ° 

K o — — * 


*0 


K. = 



K i.i 


3 q o ' £ i ? 3 £ i' £ i r 15 ( q o' e i)( c lo' e j) - 

,5 j + 2 ^5 


L % 


% 


' ~~2 


9o 


K i.i. k = 


o e • e . 

J i j 

2 <4 


15 <3o-ei)(3o*«j)^ 15 3 0 -«i «j -6 k _ 35 (5o-«i)(q 0 -ej)(qo- ? k) 


<Jo 


<lo 


+ T 




q o 


K. . . , 

i ,j ,k ,1 


is e i ' £ j e k • 7 i 105 < q o* ;Ji 315 (q 0 ‘ £ i>(5o‘ «j)(qo' ? k)(3<r £ ;) 

8 a 7 + 4 ~1T a ii 

q 0 Qq q 0 


15 (q 0 - £ i)(«j • « k ) _ 35 ( 5 o- ? i)(q 0 - ^XSo' s) 

~2 1 € , + — e 




6 I + ' 2 


<5o 


i f ’ 


where i, j , k, and l take on all positive integer values. 

This gives the expansion of -q/q 3 through the fourth order in e, and with the proper sub- 
stitutions all the terms in F through the fourth order are obtained. It is clear how these expan- 
sions are extended to any order by a straightforward process. 
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EXPANSION FOR TWO PLANETS 


In many cases the major part of the perturbations comes from the mutual attractions of the 
disturbed planet and one disturbing planet. The motions of Jupiter and Saturn are the prime ex- 
ample, and the motion of a minor planet disturbed by Jupiter is a special case. The actual com- 
putation in such cases is greatly simplified by the appearance of only two mean anomalies in the 
arguments of the trigonometric terms so that numerical double harmonic analysis can be used 
extensively. The usual series multiplications can be replaced by the more efficient process of 
evaluating the periodic terms within the double harmonic analysis routine. By holding the powers 
of the time as factors the operations are performed on functions which are periodic in two variables. 

The expansion is actually done only with respect to the secular terms, the periodic terms being 
accounted for in a direct manner. This feature provides a more complete result from the higher 
order affects of long period terms. It is necessary to make an extensive accounting for terms 
factored by rather large powers of the time. All terms through the sixth power of the time are 
given. 

Since there is only one disturbing planet in this case, in the previous section replace m. by m' , 
r. by r', and p i by p. Write 


r = r T + s T 


P = Pi + a i 


with 


Also 


Sj = Sj T + S 2 T 2 + . . . 

sj = s[ T + sj T 2 + . . . 

?i = ^ T + a 2 T 2 + . . 


5 



In this form r r , , p T , s . , s! and a. are all periodic vector functions of the two mean anoma- 

lies g and g'. The expansion of the previous section can now be used by simple substitutions. For 
example the expansion of - r/r 3 is obtained by replacing q 0 by r T , and e. by s. T 1 . In the corre- 
sponding expansions of p/p 3 and -r'/r' 3 , the terms involving the fourth power of the perturbations 
can be dropped because of the mass factor m'/(M + m ). In the solar terms 


r 




o 


put 


s = s n + s T 


^ = r 0 + s 0 , 

where s 0 contains the purely periodic terms in s. In the quantities K 0 and K i (previous section) 
resulting from the expansion of - r/r 3 in powers of Sj , it is desirable to perform a further ex- 
pansion in powers of 


S = 2- 


2 r 2 
r o 0 


The resulting quantities, being of second order, are then given with greater numerical accuracy. 
Omitting the details of the substitutions, with the disturbing force in the form 

a 2 F = a 2 f r + a 2 T + a 2 f 2 T 2 + . . . 


the first term is 


a 2 ?, =lf!i ^1 


2 r 3 
r o 


{‘-I s [‘-M 1 -! 5 - •)]} 


i/ifi-lis* fi .Is U- 

r 2 1 2 r 2 8 L 6 V 


r o + 


M+m 


Pi_ 


For the remaining terms 


r, = B~ s. . ?. * b; 7, * D* 7; * E» P, * £ [ B y > 7, * D»> * E»> S,] 

j = l 
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where 



B** 



IS a2 V 

2 r 5 
r o 


S 




11 5 

8 


B ( i)_ _3 a^S f _ 5_ s 

1 2 r* l 4 


1 




D (i) _ _ 

1 M+m _*3 

r i 


E(f) - . m ' 

1 M+m 


P 3 ! 


B ( .'~ 1) 

I 



1 


+ 3 


m J 

M+m 


r 


# 5 
I 


E ( l _1) 



BCi-2) 

i 


+ 3 


r 


5 

I 


3 

+ 2 ’ 


l 5 a 2 (r r si ) 2 


7 

I 


2 


Ch 1 00 



and 


D^" 2 ) = + 


M+m 


, a2? ;-s' 3 a 2 s ; 2 15 a*(?i-JJ>a 

O + — ■ 


2 r .s 


2 r*; 


s (i-2) _ _ m* 

1 M + m 


a 2 Pj-? 2 3 a 2 ^ 15 a 2 (Pj • (?j) 2 




+ - 


2 ,5 2 




The remaining factors D ( . lk) and k) may be deduced from b\ x ~ K) by inspection. 


(i-k> 


(i-3) _ a2? i • s 3 a 2 s, • s 2 a 2 (7 r • ijK?! • ? 2 ) 35 a 2 (7 I -i 1 ) 3 15 (7, • 1^ )s 2 

D . - + *J + o - 15 + 


2 rf 


_d-4) a2? r®4 , a2a r a 3 a 2 (7 I -7i)(?i- s 3 > 3 a 2 s 2 15 a 2 (7 t • s 2 > 2 

B = + 3 + 3 - 15 + _ ■ 


2 r 5 2 r 7 
r i r i 


105 


a ( r i‘ s l) C r I * S 2 5 15 ( r i ’ s 2 ) S 1 ( r I * S l)( S l * S 2 ^ 


2 ^ 


- 15 


w (i-5) _ - a2? i * ^5 . a2? i ■ s 4 aVj * s l )(7 I • s 4 ) 

-D ■ — + O + O — 15 


a 2 V ® 3 , e a 2 ( ? r®2)< ? r® 3 ) 105 a 2 ( ? i - ®i) 2 ( 7 r a 3 ) 105 a2 (V 7 i)( 7 r * 2 ) 2 


+ 3- 


15 


15 ( 7 I ' 7 i ) s 2 ( 7 i • s 2 )(s 1 • Sj) 15 ( 7 !-s 3 )s 2 (fj • S,)(Sj • S 3 ) 


2 rj 


- 15 


2 r; 


- 15 


= 0 , 


' a r • s. 

D* =— • 3 £ 1 

1 M+m r i5 


F * _ ■■■ 

1 " "M+m 


a '°r g i 

p\ 


n * 3 a2s i is a *( 7 r ®i ) 2 
2 = + 2 


r 5 2 r 7 
r I r i 


n' [ 3 a2s l 2 I 5 a2 (rj's;) 2 


M+m 1 ' 5 2 _/5 2 '7 

l r i r i r i 
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_* ra' f aVr? 2 3 a2a i 15 a2 (?r ? i) 2 ' 

L 9 7» - T- r * 


M+m 


P\ 


2 p| 2 


p\ 


a " s r s 2 1C a2? r s i r i • s 2 i 5 a 2 ? r s i s i 35 a 2 ( 7 r ¥ i) 


b: = + 3 




15 


2 r; 


2 r; 


and 


D *=^i 3 - 


2 F;-s' a 2 ?; -s' a^-s[?J-5' 15 a 2 7J • sjsj 2 35 a 2 (?J • sj) 2 ' 


i* * 5 


+ 3 


r'5 


15 


2 rj 7 


The remaining factors, D* and E*, can be supplied by inspection from b* taking care not to omit 
the terms factored by s! and . 


= 3 - 


- 15 


* 2 (Ji ■ s^C?! * S 3 ) 3 a 2 s 2 


2_5 2 r 7 

r I r J 


a 2 7 1 -s 1 s -s 2 15 a ! Vv| 105 a2 ( ? r a i) 2 7 r a 2 

15 - + — 


2 r? 


is a S I ; 105 a ( 7 r s i) s ? 3 i 5 a ( r r s i> 

8 r 7 + 4 r 9 

r i r i 


8 r ii 

r i 


Q * _ a S l* S 4 1C 3 ( r i ‘ S !)( r i * S 4> _ 3 S 2 ■ S 3 1C 3 ( r i ’ S 2 )( r i * S 3) 

JJ r — J ^ Jl J - ■ ■ ■ . “|" O 1 J 

5 7 5 7 

r i F I r i r l 


- 15 - 


2 r .c \2, 


a ‘ r r s i s r s 3 i5 ar r s 3 s i io5 a < r r s i) r r s 


2 vl 


2 r ? 


, - a 2 ? i , ®a»r a 2 i 5 a 27 r®i s l 105 a2 ( 7 r a 2 ) 27 r 

15 - — + 

t-7 2-7 2 -9 


15 a2s ?®i • a 2 105 a2 ( 7 i • a i) 2? i • a 2 105 a2 ( 7 i • ®i)( 7 i • ®2> s ? 
2 .7 + 2 ,9 + 2 r 9 


315 a2 ( 7 i‘ a i) 37 r a 2 
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Finally, 


a 2 (rj * s 1 )(r I * s $ ) ^ a 2 i * 2 • s 4 


4 

4 

4 


1t . a2 (4 • 4><4 • 

4) 3 a2 4 15 

a 2 (4 • s 3 ) 2 


4 

1 2 r 5 2 

r i 

4 


a 2 ?', • s, • s 4 

is a2 4 • 44 

105 a 2 (4 ' 4> 2 4 ' 4 


' 

2 4 

2 4 


a 2 7 • s s • s 
_ i<; 1 12 J 

a 2 r T • s,s, * s 

- 15 

2 . IS a2?I ■ V * 2 t 

a2 (4 • 4) 4 ■ 4 4 • 4 

4 

4 

4 

4 

15 a2 4 ' 4 4 

35 a2 ( v r 44 



2 4 * 

2 4 



15 a 2 s 2 s 1 -5 3 

105 a 2 (V4> 2 4- 

4 105 a2 4 - 44 ’4 s 2 

315 a 2 (4’4> 3 4’4 

" 2 4 * 

2 4 

4 

2 - r ii 

1 I 


15 

4 


a2 44 


15 a 44 ' 44 105 a 44 ~ 44 s 2 + 105 a 44 ' 44 s 2 

2 7 + 4 9 4 9 

r I r I 1 


4 - 105- 


945 a 44 ■ 44 (4 ' 44 
4 .11 


In forming the D and E factors, the terms involving the fourth power of the perturbations may be 
dropped since the mass factor will reduce these to the fifth order. 

These are the final formulas in the form in which they have been programmed except that the 
fourth order terms have not yet been included in a 2 f 5 and a 2 f 6 . The importance of the higher 
powers of the secular terms depends on the time interval for which the results are to be accurate. 
The fourth and higher order terms make no contribution to f r , or f 2 and the fifth order terms 
will make no contribution to f 3 so these formulas are complete to all orders. If fifth order terms 
are added then f 4 will be complete, etc. In these considerations the order is given by the power 
of the secular terms which appears, the mass factor m '/(M + m) increasing the order by one. The 
essential difference from the usual approach is that here all of the terms in the perturbations are 
taken to be of the first order. It is only assumed that the terms for higher powers of the time 
decrease so that the process can be terminated at some point, and this assumption is satisfied in 
the applications with an appropriate time scale. 
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The expansion was made as if the terms s 0 , s p s 2 , , , . in the perturbations were already 
known. A method of successive approximations, or iteration, can be used starting with any approxi- 
mate solution. If at first s and s' are put equal to zero then all terms in the expansion are zero 
except f x and 


a 2 F - a2m ' 'll _Ii_ 

M+ m n 3 J3 

Jo o _ 

which is the usual first order expression. Upon integrating, approximate values are obtained for 
s 0 and Sj (Reference 2). If Musen's method is used, the perturbations can be obtained in the form 


Then s 0 is obtained from the periodic terms in <*■ , P and y , and is obtained from the secular 
terms etc. If (3 contains large long period terms, it is at this point that the constant term and the 
term with twice the long period argument in a are adjusted (Reference 5). Letting 0 be the 
long period argument, f3 contains the terms 

C cos 6 + S s i n 0 , 


and - f3 2 /2 contains the terms 


-i (C 2 + S 2 ) (C 2 - S 2 ) cos 2 e -1 CS sin 2 6. 

When the eccentricity is neglected, these terms give the main differences between a and Hansen’s 
component v . If these terms from the first order are added to a, the results obtained in the next 
approximation are much improved, as found from experience. Furthermore, in the higher approxi- 
mations it is worthwhile to make the corrections corresponding to the changes in the coefficients 
C and s. In the case of Jupiter and Saturn this process will reduce the number of necessary 
iterations by about one half. In any particular case the importance of these corrections depends 
on the size of the long period terms. 

The formulas given in this section have been found to be the most convenient for computing 
general planetary perturbations in rectangular coordinates. Several others have been tried. All 
the terms in the mutual attractions of two planets are computed using a single basic program with 
separate subroutines for evaluating the factors of the various powers of the time. The smaller 
cross- action terms which are functions of three mean anomalies are discussed in the next section. 
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CROSS-ACTION TERMS 


For the mutual attractions of two planets an extensive development of the disturbing force, F, 
has been given. Those formulas are sufficient for computing the separate affects of the individual 
disturbing planets. What remains is to give the cross actions of the separate contributions. In 
this category, for example, would be the modifications of the perturbations of Jupiter due to Saturn 
arising from the perturbations of Jupiter and Saturn due to Uranus. In addition, the solar attraction 
on Jupiter is modified by the interaction of the perturbations of Jupiter due to Saturn with the per- 
turbations of Jupiter due to Uranus. In the higher orders the verbal description of the affects 
becomes more involved, but the formulas can be derived to any prescribed degree of accuracy by 
a straightforward and well defined procedure. Either Musen T s approach (Reference 3) using 
symbolic operators or the binomial expansion method may be used. The point here is not to enter 
into an extensive expansion, but rather to give a form which has been found most suitable for 
computation. 

In the mutual attractions of two planets complete advantage was taken of the powerful method 
of numerical double harmonic analysis by evaluating the trigonometric terms within the routine. 
Considering the periodic functions of three variables which occur in the cross action terms, one 
is not faced with any theoretical difficulty in performing numerical triple harmonic analyses, but 
rather the difficulty is that the operations cannot be performed within the computer memory. There 
are too many terms. The job could be done in parts but then there is no apparent advantage over 
the use of series multiplications. 

The method given here takes a third approach wherein full use is made of double harmonic 
analysis, and the terms containing the third argument are held as factors until the final step. 

It will be sufficient at first to give only those cross action terms which are of the second order. 
The term "disturbed planet" refers to that planet for which the disturbing force is being given. 

The notations are, for the present purpose, 

r 0 the position vector of the disturbed planet as obtained from the elliptic motion 
s 0 the perturbations of the disturbed planet 
m Q the mass of the disturbed planet 

r . the position vector of the i th disturbing planet as obtained from its elliptic elements 
Sj the perturbations of the i th planet 
m. the mass of the i th planet 

p. the position vector of the i th planet with respect to the disturbed planet, p. = r. - r 0 
M the mass of the Sun. 
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With this notation the disturbed planet is also the zero th planet, 
in F are 


The second order terms 


F = 3 




15 <?o ' ®o > 2 - 

2 r 7 F| 
r o 


+ 


5 


m. 

i 

M+m 0 




1 

p i + ~ 
Pi 


p i ’ S i - 
3 l p i 
P\ 





where N is the number of disturbing planets. This formula gives all the terms of the second order, 
whereas only those involving cross actions are needed here. Put 


N 



j = 0 


j/ i 


i = 0 , 1 , 2 N 


where s . represents the perturbations of the i th planet due to the j th . 

Now substitute these expressions for the perturbations into the formula for F . 





+ 










+ 



"O.j 



15 (Vso.iKiWo.j) 

2 ~7 

r o 



+ 


N 

L 

i J=1 


m. 

i 

M+m 0 







The terms in the single summations are taken into account in the previous section, and the double 
summations give the second order cross-action terms which are to be considered. These latter 
terms are most conveniently separated into classes according to the perturbations by which they 
are factored. The total is given by 


+ F n + F m 
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where 


F, = 




o.j 


and 


N 



i.j = l 
i/j 


N 



i ,j = 1 
i^j 


the dyadics being given by 


Hj‘> = 3 r °' 50,1 


j , S 0 ,i r 0 , r o s o.i 

1+3 + 3 


lc r o* s 0,i - - 
15 r o r o , 


l (i) = ^ J_ J_ I + 3^1 

11 " M +% | p s 


and 


H 


(i) 

III 





the unit dyadic being denoted by I . The practical value of these forms comes from the fact that 
the dyadics, H, are independent of the mean anomaly of the j th planet. It is best to carry the 
development further and consider the components of F as they enter the formulas for the pertur- 
bations in Musen T s method (Reference 2). In dimensionless form with time as the independent 
variable 


and 


0 = 


y = 


i 

i 

i 


(MjD + M 2 E)d(nt) , 


(M 4 - 2a)d(nt) , 


M 3 Dd (nt) , 


14 


rOl 


to which the contributions from the constants of integration are to be added. The integrands 
obtained from the formulas: 


M x = S t • (a 2 F) , 


M 2 = S 2 • (a 2 F) , 


M, = — R • (a 2 F) , 


■J 


M. = M d(nt) , 


- P cos e 


+ Q 


sine 

/TT? 


S 2 


- P sin e + Q 


cos e - e 

✓ l - e 2 


and 


D = — [sin (rj - e) - 
L 


e sin t) + e s 


i n ej 


E = 2 — [l - cos (r) - e)] , 
r o 


where the notations are 

a the semi -major axis of the reference ellipse 
e the eccentricity of the reference ellipse 
n the mean motion in the reference ellipse 
e the eccentric anomaly in the reference ellipse 

V ' the eccentric anomaly which is held constant during the integration and replaced by e 
afterward 

, R the standard unit vectors of the orbit. 


The perturbations are given by 


s o “ a o r o + w + ^0 a ® » 


where 


~ _ 1 d7 o 
n dt 


For purposes of computation , 


r Q = aP (cos e - e) + a /l - e 2 Q sin e 


:_L [-aPsine+a /l - e 2 Q cos ej 


r Q □ a (1 - e cos e) . 

With these formulas the contributions of a 2 Fj , a 2 F n ^and a 2 F m to the components M a , M 2 and 
M 3 can be given in a convenient form. Only and M 2 are given since M 3 , which is of lesser im- 
portance, can be handled in the same way. 

Class I Terms 

The perturbations, s Q } will be in the form 


>0.J = a o,i 7 o + 00, i 5 + y o,i aS • 

To hold the multiplications by terms involving g. until the last, put 

^1 = ^1 ,1 a 0 ,j + ^1 , 2 ^0 , j + ^1 .3 ^0 , j 

and 

^2 = ^2,1 a 0,j + ^2,2 ^0,j + ^2,3 ^O.j ' 

The factors are given by the formulas 


^1 ,1 " * r o » 


Kj 2 qLj* w 
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1 


Kj 3 - Lj * aR , 
K 2 ,i “ ^2 ‘ r o ’ 


and 


with 


K 2i2 =L 2 - w , 

K 2 3 = L 2 * aR , 


a 2 Hj i} = a 2 



and 


L. = S • a 2 H 


(O 





Class II Terms 

Terms in this class are computed from the same formulas as the terms in Class I except that 


and 


L i ~ s i 


i (i) 
1 ri 




Class III Terms 

The perturbations of the disturbing planets will be in the form 

s. . = a. . 7 . + 3. . w. + y. .a. R. 

l ,j l.j l 'l ,} i 'l.J i l 

corresponding to those of the disturbed planet. Now put 

^1 “ ^1 ,1 tt i ,j + K 1 ,2 A ,j + ^1 .3 ,j 1 
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M 2 = K 2,i a i.i + K 2.2A.i + K 2.3^i,j. 


K i.i =L r 


V 1 ,2 


V w i 


and 


with 


k i , 3 = L r a ; R i • 


K 2,l = V r i ■ 


K 2,2 *>V W i • 


*2.3 = V a i R i 


Lj = S, 


1 (1) - 
1 III “ 


M + m 0 


S, -3- 


■ p. + 3 ■ 


Pi 


and 


a 2 


M + m n 



p. + 3 


In each class of terms the factors Kg m depend only on the two mean anomalies g 0 and g., so 
these factors can be expanded using double harmonic analysis. After these factors are formed, 
the coefficients in a 0 . , /3 Q . , y 0 . or a £ > . , , y. . are taken in groups corresponding to fixed 

multiples of the mean anomaly g. . The resulting terms in M 1 and m 2 are obtained by series 
multiplications, and the terms in a and /3 are computed by the integration procedure and written 
out. Then the next group of terms is read into the computer. Using this procedure it is possible to 
work entirely within the computer memory, and this makes the execution much more efficient than 
would be possible otherwise. In some cases it may be desirable to combine the terms which have 
been separated into classes I and II, and this can easily be done. 

The discussion has been limited to terms of the second order, but it should be clear how the 
same method can be applied to many of the most important higher order terms by modifying the 
vectors L 1 and L 2 to include powers of s 0 . and s i 0 and by taking terms with three or more 
arguments in the factors. 


18 



ACKNOWLEDGMENT 


The author expresses his appreciation to Dr. Peter Musen who suggested the work on general 
perturbations, to Neal Paris for his assistance in programming, and to Edson Goodrich for checking 
the formulas. 

(Manuscript received May 3, 1965) 


REFERENCES 

la. Hill, G. W., "A New Theory of Jupiter and Saturn,” Astronomical Papers Prepared for the 
Use of the American Ephemeris and Nautical Almanac, 4:1-576, 1890. 

lb. Hill, G. W., "A New Theory of Jupiter and Saturn," Memoir No. 50 in: "The Collected 
Mathematical Works of George William Hill, Volume 3," Washington: Carnegie Institute of 
Washington, Publication No. 9, January 1906. 

2. Musen, P., and Carpenter, L., "On the General Planetary Perturbations in Rectangular 
Coordinates,” J. Geopliys. Res. 68(9):2727-2734, May 1, 1963. 

3. Musen, P., "On the General Perturbations of the Position Vectors of a Planetary System," 
NASA Technical Note D-2692, 1965. 

4. Carpenter, L. H., "Computation of General Planetary Perturbations, Part I," NASA Technical 
Note D-1898, August 1963. 

5. Carpenter, L. H., "Computation of General Planetary Perturbations, Part II. A Comparison 
of Components," NASA Technical Note D-2852, June 1965. 


NASA- Langley, 1966 


19 


"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless 

of importance as a contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: Information receiving limited distri- 

bution because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in con- 

nection with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 

language considered to merit NASA distribution in English. 

TECHNICAL REPRINTS: Information derived from NASA activities 

and initially published in the form of journal articles. 

SPECIAL PUBLICATIONS: Information derived from or of value to 

NASA activities but not necessarily reporting the results of individual 
NASA-programmed scientific efforts. Publications include conference 
proceedings, monographs, data compilations, handbooks, sourcebooks, 
and special bibliographies. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


